Theoretical description of the double beta decay of 160 Gd 
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The half-life of the (3(h v decay of 160 Gd, a process which was previously reported as theoretically 
forbidden in the context of the pseudo SU(3) model, is estimated by including the pairing interaction. 
Thus, different occupations can be mixed by the interaction opening new channels for the decay. 
Explicit expressions are presented for the mixing induced by the pairing force. Matrix elements for 
the Pfav and (3(3o v decays are calculated, in the present pseudo SU(3) approach, by assuming both 



a dominant component in the wave function of the ground state of Gd and a more general model 
space. The results, of the calculated (5(3 half-lives, suggest that the planned experiments would 
succeed in detecting the (3(32i> decay of 160 Gd and, eventually, would improve the limits for the zero 
neutrino mode. 
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The flux of solar and atmospheric neutrinos has been measured with increasing precision, and the data offer direct 
evidence on the presence of neutrino oscillations [1, 2]. These findings imply that at least neutrinos of some flavor 
should be massive. The difference between the square of the masses of neutrinos belonging to different families can 
be extracted from the experimental data, relying on theoretical models of mass hierarchies and textures [3]. On the 
other hand, their absolute scale cannot be obtained from these experiments. 

The neutrinoless double beta decay (P0ol>)i if detected, would provide the complementary information needed to 
determine neutrino masses, and would also offer definitive evidence that the neutrino is a Majorana particle, i.e. that 
it is its own antiparticle [4, 5]. 

Theoretical nuclear matrix elements are needed to translate experimental half-life limits, which are available for 
many /3/3-unstable isotopes [7, 8], into constrains for the effective Majorana mass of the neutrino and, eventually, 
for the contribution of right-handed currents to the weak interactions. Thus, these matrix elements are essential to 
^ ' understand the underlying physics. 

The two neutrino mode of the double beta decay {fifiiv) is allowed as a second order process in the standard model. 
It has been detected in ten nuclei [7, 8] and it has served as a test of a variety of nuclear models [5]. The calculation of 
the /3/?2i/ and (3(3§ v matrix elements requires the use of different theoretical methods. Therefore a successful prediction 
of the former cannot be considered a rigorous test of the latter [6]. However, in most cases it is the best available 
proof we can impose to a nuclear model used to predict the (3(3q u matrix elements. When possible, the test should 
also include the calculation of the energy spectra, electromagnetic transitions and particle transfer observables, in the 
neighborhood of the double beta emitters. 

Many experimental groups have reported measurements of (3(3 processes [7, 9]. In direct-counting experiments the 
analysis of the sum-energy spectrum of the emitted electrons identifies the different /3/3-decay modes [10] . 

The pseudo SU(3) approach has been used to describe many low-lying rotational bands, as well as BE(2) and B(M1) 
intensities in rare earth and actinide nuclei, both with even-even and odd-mass numbers [11, 12, 13, 14, 15, 16]. The 
(3(3 half-lives of some of these parent nuclei to the ground and excited states of the daughter were evaluated for the 
two and zero neutrino emitting modes [17, 18, 19, 20, 21] using the pseudo SU(3) scheme. The predictions were found 
to be in good agreement with the available experimental data for 150 Nd and 238 U. 

Based on the selection rules of the simplest pseudo SU(3) model, the theory predicts the complete suppression of 
the (3(3 decay for the following five nuclei: 154 Sm, 160 Gd, 176 Yb, 232 Th and 244 Pu [17]. It was expected that these 
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forbidden decays would have, in the best case, matrix elements that would be no greater than 20% of the allowed 
ones, resulting in the increase, by at least one order of magnitude, of the predicted half-life [21]. Experimental limits 
for the (3(3 decay of 160 Gd have been reported [22, 23]. Recently it was argued that the strong cancellation of the 2v 
mode in the (3(3 decay of 160 Gd would suppress the background for the detection of the Ov mode [24]. 

In the present contribution we extend the previous research [17, 18, 19, 20, 21] and evaluate the (3(3 half-lives 
of 160 Gd using the pseudo SU(3) model. While the 2v mode is forbidden when the most probable occupations arc 
considered, states with different occupation numbers can be activated by the pairing interaction. The amount of 
this mixing is evaluated, and the possibility of observing the (3(3 decay in 160 Gd is discussed for both the 2v and Of 
modes. The analysis is performed firstly for the extreme case of a single active configuration in the initial nucleus, 
and afterwards, the results are compared with the ones obtained by enlarging the model space. 

The paper is organized as follows. In Section II the pseudo SU(3) formalism and the model Hamiltonian are briefly 
reviewed. In Section III the 160 Gd and 160 Dy ground state wave functions are built. Section IV and V contain the 
explicit formulae needed to evaluate, using the pseudo SU(3) scheme, the (3(32 V and (3(3q v matrix elements, respectively. 
In Section VI the 160 Gd (3(3 nuclear matrix elements and half-lives are presented. Conclusions are drawn in the last 
Section. 



II. THE PSEUDO SU(3) FORMALISM 

In order to obtain a microscopic description of the ground states of 160 Gd and 160 Dy we will use the pseudo SU(3) 
model, which successfully describes collective excitations in rare earth nuclei and actinides [11, 12, 13, 14, 15, 16, 25], 
and it has been used to evaluate the half-lives of the (3(32 V decay to the ground and excited states and (3(3q v decays 
of six heavy deformed nuclei [17, 18, 19, 20, 21], and the double electron capture decays in other three nuclei [26]. 

In the pseudo SU(3) shell model coupling scheme [27] normal parity orbitals with quantum numbers {f],l,j) are 
mapped to orbitals of another harmonic oscillator with (i) = rj — 1,1, j). This set of orbitals, with j = j = I + s, 
pseudo spin s = 1/2 and pseudo orbital angular momentum I, define the so-called pseudo space. For configurations 
of identical particles occupying a single j orbital of abnormal parity, a simple characterization of states is made by 
means of the seniority coupling scheme. 

The first step in the pseudo SU(3) description of a given nucleus consists in finding the occupation numbers for 
protons and neutrons in the normal and abnormal parity states n% ,n^,n^ [17]. These numbers are determined 
filling the Nilsson levels from below, as discussed in [17]. 

For even-even heavy nuclei it has been shown that if the residual neutron-proton interaction is of the quadrupolc 
type, regardless of the interaction in the proton and neutron spaces, the most important normal parity configurations 
are those with highest spatial symmetry {f a } — {2™°/ 2 } [25]. This statement is valid for yrast states below the 
backbending region. It implies that = S v = 0, i.e. only pseudo spin zero configurations are taken into account. 

Additionally in the abnormal parity space only seniority zero configurations are taken into account. This simplifi- 
cation implies that J A = J(} = 0. This is a very strong assumption quite useful in order to simplify the calculations. 

Many-particle states of n a active nucleons in a given normal parity shell i] a , a — v or n, can be classified by the 
following chains of groups: 

I 1 ™"} {fa} {fa} la (A a ,Ma) $a K a L a 

U(Sl% ) D U(Sl% /2) x U(2) D SU(3) x SU(2) D SO(3) x SU(2) D SUj{2), (1) 

where above each group the quantum numbers that characterize its irreps are given and j a and K a are multiplicity 
labels of the indicated reductions. The most important configurations are those with the highest spatial symmetry 
[25, 28], namely those with pseudo spin zero. 

The model Hamiltonian contains spherical Nilsson single-particle terms for protons (H sPi „) and neutrons (H sp ^), 
the quadrupole-quadrupole (Q ■ Q) and pairing interactions (V pa i r ), as well as three 'rotor-like' terms which are 
diagonal in the SU(3) basis. 

H = H sPy7T + H sPtV + Vpair — ^ X Q ' Q + a Kj + b J 2 + A asym C^.- (2) 

This Hamiltonian can be separated into two parts: the first one includes Nilsson single-particle energies and the 
pairing and quadrupole-quadrupole interactions (Q is the quadrupole operator in the pseudo SU(3) space). They 
are the basic components of any realistic Hamiltonian [29, 30] and have been widely studied in the nuclear physics 
literature, allowing their respective strengths to be fixed by systematics [29, 30]. In the second one there are three 
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rotor terms used to fine tune the moment of inertia and the position of the different K bands. The SU(3) mixing is 
due to the single-particle and pairing terms. 

The three 'rotor-like' terms have been studied in detail in previous papers where the pseudo SU(3) symmetry was 
used as a dynamical symmetry [25] . In recent works, a, b, and A asym were the only parameters used to fit the spectra 
[13, 15, 16, 31, 32]. 

The spherical single-particle Nilsson Hamiltonian is 



H sp = Hlj (t]+ -) - Khoj {2l ■ s + M^} = ^2 e '3i) a \ a i ( 3 ) 



with parameters [29] 



hujo = 41A _1/3 [MeV], = 0.0637, k v = 0.0637, (4) 
Mtt = 0.60, fi u = 0.42, 

The pairing interaction is 



Vpair = ~ \G ^ a j a ] a J' Q j' ( 5 ) 



3,j' 



where j denotes the time reversal partner of the single-particle state j, and G is the strength of the pairing force. 
In principle we can start with an isospin invariant pairing force, which contains proton-proton, neutron-neutron and 
proton-neutron terms with equal strengths. However, in the limited Hilbert space employed there are not protons and 
neutrons in the same orbitals, and the proton-neutron pairing term is ineffective. In recent works [13, 15, 16, 31, 32], 
the pairing coefficients G^m were fixed following [29, 30], with values: 

G n = ^- = 0.132 MeV, G v = ^ = 0.106 MeV. (6) 



III. THE GROUND STATE OF GD AND 160 DY 



In this section we shall present the relevant results of the pseudo SU(3) model, for the wave functions of the 
participant nuclei. We shall assume, in subsection (III. A), that the wave function describing the ground state of 
160 Gd has only one component and that the corresponding wave function for 160 Dy has two components connected 
by the pairing interaction. This assumption allows us to perform a simple estimation of the double beta decay, as 
shown in sections IV and V. In subsection (III.B) we shall include other configurations in the model space, to asses 
the quality of the simplified description. 



A. The one-component model 



With a deformation of e = 0.26 [33], the most probable occupations for the 160 Gd 14 valence protons are 8 particles 
in normal and 6 particles in unique parity orbitals, and for the 14 valence neutrons are 8 particles in normal and 6 
particles in unique parity orbitals. 

After the detailed study of 156 . 158 i 1 60q c j i so t p es performed in [16, 32], it becomes clear that the pseudo SU(3) 
model is a powerful tool in the description of heavy deformed nuclei. Up to four rotational bands, with their intra- 
and inter- band B(E2) transition strengths, as well as B(M1) transitions, were reproduced, with a very good general 
agreement with the experiment. 

The dominant component of the 160 Gd ground state wave function [16] is 

| 160 Gd,0+) = | (h u/ 2)l, = 0; (*i3/ 2 ) 6 „ = 0)a (7) 

|{2 4 } T (10,4) ;r ;{2 4 } I ,(18,4) ly ; 1(28, 8)K = 1, J = 0) N . 

In the first series of papers [17, 18, 19, 20, 21] we have evaluated the matrix elements by taking into account only 
the leading SU(3) coupled proton-neutron irrep. In recent publications it was shown that the leading irrep represents 
approximately 60 % of the wave function in even-even Dy and Er isotopes [16, 31]. 
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Assuming a slightly larger deformation for 160 Dy, the most probable occupations for 16 valence protons are 10 
particles in normal and 6 particles in unique parity orbitals, and for the 12 valence neutrons are 6 particles in normal 
and 6 particles in unique parity orbitals. A detailed study of 160 > 162 > le4 Dy i so t pes has been performed in [32]. 

The dominant component of the wave function is 



|160 



Dy,0+(a)> 



| (h u /2)l, J? = 0; (* 13 /2)«, # = 0}a 
|{2 5 K(10,4),;{2 3 },(18,0) I ,; 1(28, 4)K = 1, J = 0) N . 



The two neutrino double beta operator annihilates two neutrons and creates two protons with the same quantum 
numbers r], I. It cannot connect the states | 160 Gd,0 + ) and | 160 Dy, + (a)) and the transition becomes absolutely 
forbidden. This is the selection rule found in [17]. 

However, the pairing interaction allows for the mixing between different occupations. In the deformed single particle 
Nilsson scheme it takes an energy AE to promote a pair of protons from the last occupied normal parity orbital to 
the next intruder orbital. This excited state has 8 protons in normal and another 8 protons in unique parity orbitals, 
and its wave function has the form 



1 160 



Dy,0+(6)> 



| (h u/2 )l, J? = 0; (i 13/2 )«, J? 
|{2 4 } 7r (10, 4),; {2 3 },(18, 0),; 1(28, 4)K = 



-0) A 
1,J = 



0) 



N- 



The two neutrino double beta decay of 160 Gd can proceed to this state. 

As a first approximation, we shall describe the 160 Dy ground state as a linear combination of these two states: 



|160t 



Dy, 0+) = a | ibu Dy, 0+(a)> + b | lbu Dy, 0+(6)} 



(8) 



with |a| 2 + |6| 2 = 1. 

The only term in the Hamiltonian (2) which can connect states with different occupation numbers in the normal 
and unique parity sectors is the pairing interaction. In the present case, the Hamiltonian matrix has the simple form 



H 







L pa%r 



IL pair 

AE 



(9) 



with 



b pa%r 



( 160 Dy,0+(6)|T/ pair | 160 Dy,0+(a)) 



(-1)^ +1 G, 



yj « + 2)(2# + l-n£)Y,\/ + 1) 



J2 ((0^)l[ T ,(0r? 7r )l/ 7r ||(A O Mo)10) 

(Aojuo) 

(A-,<) (A ,Mo) P, 

{K^D (o,o) (xlfil) i 

{\ a ,^) (A ,Mo) (A 6 /x 6 )cr p 
1 1 1 



^((A a , M a )l 0,(A o /io)l 0||(AV b ).10) P 



III [ao^^o^iP '^ III M,ti)) P „ 



(10) 



In the above formula (..,..||,,) denotes the SU(3) Clebsch-Gordan coefficients [34], the symbol [...] represents a 
9 — A/i rccoupling coefficient [35], and (.. ||| .. ||| ..) is the triple reduced matrix elements [18]. 
The lowest eigenstate has an energy 



E 



AE 



i - \ i 



(2h 



pair 



V AE 



and the components of the 160 Dy ground state wave function are 



E 



(11) 



(12) 



It is a limitation of the present model that AE has to be estimated from the deformed Nilsson single particle mean 
field, instead of the evaluation of the diagonal matrix element of the Hamiltonian (2). The use of seniority zero states 
for the description of nucleons in intruder orbits inhibits the direct comparison of states with different occupation 
numbers. A formalism able to describe nucleons in both normal and unique parity orbitals in the same footing is 
under development [36]. The effect AE has on the double-beta-decay half- lives is studied in Section VI. 
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B. The enlarged model space 

In addition to the pseudo SU(3) irreps described in the previous subsection, one has, in an enlarged model space, 
other potential components to the ground state wave functions of 160 Gd and 160 Dy. In Tables I and II we are listing 
the four configurations included in the description of 160 Gd and the ones for 160 Dy, respectively. 



state 




ni 


JV 

n v 




(Ajt , /i-TT ) 


(A„,//„) 


(A,M) 


C(A, M ) 


ii 


8 


6 


8 


6 


(10, 4) 


(18,4) 


(28,8) 


0.7253 


12 


8 


6 


10 


4 


(10, 4) 


(20,4) 


(30,8) 


0.4848 




10 


4 


8 


6 


(10, 4) 


(18,4) 


(28,8) 


0.4066 


u 


10 


4 


10 


4 


(10, 4) 


(20,4) 


(30,8) 


0.2713 



TABLE I: The four configurations included in the description of Gd. The number of proton s and neutrons, in normal and 
in intruder orbits and the corresponding pseudo SU(3) irreps are listed. The last column shows their amplitudes, in the wave 
function of the ground state of 160 Gd. 

The wave function, of the ground state of 160 Gd, is written 

k 

and the one of 160 Dy, is given by 

1 160 Dy) = Y,Ci S) \f k ), 

k 

where the amplitudes and Cjf^ are listed in the last column of Tables I and II, 



state 








n$ 




(A„, fi v ) 


(A,M) 


C(A, M ) 


fi 


10 


6 


8 


4 


(10, 4) 


(18,4) 


(28,8) 


0.7639 


h 


10 


6 


6 


6 


(10, 4) 


(18,0) 


(28,4) 


0.5446 


h 


8 


8 


8 


4 


(10, 4) 


(18,4) 


(28,8) 


0.2690 


U 


8 


8 


6 


6 


(10, 4) 


(18,0) 


(28,4) 


0.2181 



TABLE II: The four configurations included in the description of Dy. The configurations and amplitudes are listed following 
the notation given in the captions to Table I. 

(9), for the configurations of Tables I and II, is a 4x4 matrix, with matrix elements h pa i r of the form (10). The 
diagonal matrix elements are the energies needed to promote a pair of protons, neutrons or both from a normal to 
an intruder parity orbital or viceversa. They are estimated from the deformed single particle Nilsson diagrams. The 
values of AE, for the case of 160 Gd are 0.54 MeV, 0.81 MeV and 1.35 MeV. The corresponding quantities for 160 Dy are 
of the order of 0.54 MeV, 1.71 MeV and 2.25 MeV. The diagonalization of the mixing matrix yields the amplitudes 
C of the wave functions, and they are shown in the last column of Tables I and II. As one can see from these tables, 
the mixing induced by the pairing interaction is sizeable although the dominant component of each wave function is 
still the irrep considered in the restricted model space. 

Naturally, and in order to estimate the effect of this mixing upon the double beta decay process, we should 
compute explicitly the corresponding matrix elements, which are the sum of products of amplitudes and individual 
matrix elements. It will be done in the following sections. 

IV. THE /3f3 2v HALF-LIFE 

The inverse half- life of the two neutrino mode of the /3/3-decay, (302v, can be cast in the form [37] 

1 = G 2v | M 2u | 2 . (15) 
where G 2v is a kinematical factor which depends on Qpp, the total kinetic energy released in the decay. 



(13) 
(14) 

respectively. The mixing matrix 



7# a (0+-0+) 
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The nuclear matrix element is written 

rT (o+ ||r|| 1+) (i+ ||r|| o+) 

M 2v « M£ T - V W " , (16) 

with the Gamow- Teller operator V expressed as 

r m = a mst~ = J2 a ^'^K^b^ "v*l* hiiJ™> m=l,0,-l. (17) 

The energy denominator is /jjv = Ej + En — Ei and it contains the intermediate En, initial Ei and final Ef energies. 
The kets |1^) denote intermediate states. 

The mathematical expressions needed to evaluate the nuclear matrix elements of the allowed g.s. — ► g.s. (3(3 decay 
in the pseudo SU(3) model were developed in [17]. Using the summation method described in [17, 38], exploiting 
the fact that the two body terms of the SU(3) Hamiltonian commute with the Gamow- Teller operator (17) [18], 
resumming the infinite series and recoupling the Gamow- Teller operators, the following expression was found: 

Mg T = V3 E 2 S$S#< / I [[4 ® ~a,Y ® [al ® a,,] 1 ] J= ° | 0+) (18) 

where 7r = (77^, l w ,j v ) and f = (?7„, l v ,j v ), and _E = + m e c 2 . 

In the following we shall analyze the nuclear matrix element (18) for the (ifav decay of the ground state of 160 Gd, 
Eq. (7), to the ground state of 160 Dy, Eq.(8). Each Gamow- Teller operator (17) annihilates a proton and creates a 
neutron in the same oscillator shell and with the same orbital angular momentum. In the case of the (3(3 of 160 Gd it 
means that the operator annihilates two neutrons in the pseudo shell r\ v — 5 and creates two protons in the abnormal 
orbit h\\/2- As a consequence the only orbitals which in the model space can be connected by the (3/3 decay are those 
satisfying 77^ = r\ v = 77, that implies l n = l v = 77, j v = r] — | and j v = 77 + |. 

Under these restrictions the /3/32u decay is allowed only if the occupation numbers obey the following relationships 

n nj — + 2 , n t.j — n v,i 1 

N N N N o V ly J 

It follows that, when using the restricted configuration space in Eq. (18), only one term in the sum survives and thus 
the nuclear matrix element M<i v (18) can be written as 

MgT = ^! (0 +! [[a t M i [a t g, ^1] ^ | + >) (20) 

where the energy denominator £ is determined by demanding that the energy of the Isobaric Analog State equals the 
difference in Coulomb energies Ac- It is given by [17] 

£ = E Q + e(r7 w ,Z w , j v = j v + 1) - z{r) v ,l v ,j v ) = E - fujk^2j n + A c . 

(21) 

A c = ffi[2Z + 1 - 0.76((Z + I) 4 / 3 - Z 4 ^)]MeV. 

As it was discussed in [17], Eq. (20) has no free parameters, being the denominator (21) a well defined quantity. 
The reduction to only one term is a consequence of the restricted Hilbert proton and neutron spaces of the model. 
The initial and final ground states are strongly correlated with a very rich structure in terms of their shell model 
components. 

The low energy levels are assumed to have pseudo spin S — 0, a fact that again simplifies the evaluation of the 
above sum since L = J. 

The calculation of the matrix elements, Eq. (20), in the normal space is performed by using SU(3) Racah calculus 
to decouple the proton and neutron normal irreps, and expanding the annihilation operators in their SU(3) tensorial 
components. For the particular case of 160 Gd, in the spirit of the single model space of subsection III. A, the flfcv decay 
can reach only the second component of the wave function of the ground state of 160 Dy (8), and for this reason it is 
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proportional to the amplitude b. The expression for the matrix elements of the Pfov channel reads 



M 2 G /( 160 Gd^ 160 Dy) = b 2^1 



«+2)(2j T + l-»^) 
Jtt-1 



1/2 



E ((Ot?)1Z,(0?7)1Z||(AoMo)10>i £((28,8)1 0,(A o ,Mo)l 0||(28,4)10) p 

(A 0l( uo) P 

(10,4) (0,0) (10,4) 1 



E 

p' 



(18,4) (Ao^o) (18,0) P > 
(28,8) (AoMo) (28,4) p 
1 1 1 



((18,0) HI [ao^iao^ip^) ||| (18,4)),,. 



(22) 



In the so-called large model space, of subsection III.B, the expression for the matrix element of the f5fti v decay 
mode is 



M GT ( 160 Gd^ 160 Dy) = ECT C, 



k,l 



2j„ 



(n^ fc +2)(2^, fc + l- 



1/2 



E ((0^)n,(07?)i;||(AoMo)10)iE((Afc,Mfe)l 0,(A ,mo)1 0||(A,, W )10>„ 

(A ,/jo) P 

{K,k,^,k) (o,o) (K.hVir.i) i 

(A^/c, Mi/,fc) (AoMo) (A 
(Xk,Hk) (AoMo) (A;,m0 P 
1 1 1 



E 



(23) 



((A/,W) HI [^),i5o^),i] (AoAio) III (Afe,Mfe))p'- 



V. THE f3(3 0v HALF-LIFE 

For massive Majorana neutrinos one can perform the integration over the four-momentum of the exchanged particle 
and obtain a "neutrino potential" which for a light neutrino (m u < 10 McV) has the form 



H { r,E)=™r dq ?m 

V ' *r J q + E 



(24) 



where E is the average excitation energy of the intermediate odd-odd nucleus and the nuclear radius R has been added 
to make the neutrino potential dimensionless. In the zero neutrino case this closure approximation is well justified 
[39]. The final formula, restricted to the term proportional to the neutrino mass, is [4, 37] 



(to?)- 1 - ( 



(25) 



where Go„ is the phase space integral associated with the emission of the two electrons. The nuclear matrix elements 
M „ are [37] 



9 a 



with 



m« = <o; no* not) , 



(26) 



(27) 



where the kets |0^ ) and |0^) denote the corresponding initial and final nuclear states, the quantities gy and gA are 
the dimensionless coupling constants of the vector and axial vector nuclear currents, and 



O gt = E Ogl = E Zmt- ■ a n t-H(\f m - f n \,E) 

m.n m.n 

o F = E o F mn = E t-t-H(\f m - r n \,E) , 



(28) 
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being a the Pauli matrices related with the spin operator and t~ the isospin lowering operator, which satisfies 
t~\n) = \p). The superindex GT denotes the Gamow- Teller operator, while F indicates the Fermi operator. In the 
present work we use the effective value (^7) 2 = 1-0 [42]. 

Transforming the transition operators to the pseudo SU (3) space, we have the formal expression 



O a = 



+ 



+ 0' 



+ 



(29) 



where the subscript index NN, NA, . . . are indicating the normal or abnormal spaces of the nucleon creation and 
annihilation operators, respectively. Given that we use the Nilsson scheme to obtain the occupation numbers, we are 
considering only nucleon pairs. For this reason only the four type of transitions listed above give a non-vanishing 
contribution to the /3/?2i/ matrix elements. 

In a previous work we have restricted our analysis to six potential double beta emitters which, within the approx- 
imations of the simplest pseudo SU(3) scheme, are also decaying via the 2v [3(3 mode. They include the observed 
150 Nd — > 150 Sm and 238 U — > 238 Pu Pfov decays. In this case two neutrons belonging to a normal parity orbital decay 
in two protons belonging to an abnormal parity orbital. The transition is mediated by the operator 0\ N . Under 
the seniority zero assumption for nucleons in abnormal parity orbitals, only proton pairs coupled to J=0 are allowed 
in the ground state. The matrix elements is 



M « (A n N„) = £ Mg v {A«N v ,j v ) = - £ 



2j„+l 



( n A+2)(2j 7 , + l-ni) 



(ti£j£)J = M - U\0"\(j v j v )J = M = 0) 



2(2i„+l) V 2jV+1 
3v 3v ' v ' 

£ £ ((\l,4)K,L„, (Xl,fj,l)K' v L w \\(X f , pf)10) {{>i,ti)K„L„, (A*,, ti)K v L„\\ (A*, ^)W) 
E ((A^, (Ao, ^o)10|| (A£, Ai/)^^)^ ((0, (0, ^)1^|| (A , ^ )10) 

(A ,po)po 

)L ; i] (Ao ' M) ^ =0 IH(At,Mt)),o- 



(30) 



((a£,m£)III[>,^)L;i>^ 



We have implicitly defined M^ijy) as the contribution of each normal parity neutron state j v to the nuclear matrix 
element in the transition (j v ) 2 — > (j^) 2 ■ 

The transitions which are forbidden for the /?/?2i/ decay are allowed for the zero neutrino mode, due to presence of 
the neutrino potential, instead. In the simplest model space, the (3(3o v of 160 Gd has finite contributions for the two 
components with different occupation numbers in the 160 Dy final state. There are two terms in the (3(3q v decay: one 
to the basis state which has allowed (ifov decay, and one to the state with forbidden (Hfov decay. In the first case the 
above equation must be used. The second case involves the annihilation of two neutrons in normal parity orbitals, 
and the creation of two protons in normal parity orbitals. This transition is mediated by the operator N ^. The 
(3(3q v matrix element has the form 



Mg> v (N v N v ) 



= £ M« (N n N v , J) = -\ E V(2A + l)(2^ + l)(2> + l)(2.v + l) 

J 3^3\'3v3„' 

E V2^ + T ((jnjn>)J\O a \(j„j„>)J) W(hjy w ,J w ,j w ) W{l v j\j v ,~l v .j v ) 
J 

E ((fk, 0)14, (^,o)iiV || (xs,i4)ksj) E ((o,v„)il, (o,»?„)ii„, \\(\z,rt)K5J) 

E ((AJ,m5)^J,(A^^)^J||(A o ,mo)10) Po E((A%m010,(Ao,Mo)10||(A/, m /)10) p 

(A ,po)po P 

(Aix) (A5,mS) (H,d) p« 



(31) 



E 

PxP„ 



((4,4) 



(At,A*t) (Ag,A«) (A£,p£) P„ 
(A*, M*) (Ao,Mo) (AV) P 
1 Po 1 

((A^M^III^O,^)^;!^,^)?, 



l(AJ,pJ):S=0 



(A;,<))p„ 



1(as,mS);S=oi 



I (At, Mt)), 



In the above expression the W(..., ...) are Racah Coefficients [43]. The two-body matrix element can be expanded in 
its L,S components 



((MJ\O a \(j^)J) = 





1-K l^' 










1 1 

2 2 




hi 


2 2 ^ ' 


LS 


k Jn ]k' 






k jv jv< J t 



bb 

'22' 



bb 

'22' 



(32) 
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where the x{--} are Jahn-Hope coefficients [43] and 

T • T(GT) = ai ■ (72 T • T(F) = 1 a = GT or F. (33) 

In order to evaluate the spatial matrix elements, we have used the Bessel-Fourier expansion of the potential [44], 
which gives 

{{hh)LM\H{r)\{hh)LM) = Y J {-l) h+u+L (21 + 1) (h\\Ci\\h) (i 2 ||Ci||i 4 ) WQihhk; LI) R\hh,hh) (34) 

where (li\\Ci\\lj) arc the reduced matrix elements of the unnormalizcd spherical harmonics C; m (f2) = ^^Yi m (fl) 

and R l (lil2, hh) are the radial integrals described in Appendix A of Ref. [18]. They include the effects of the finite 
nucleon size and the short range correlations, as explained in Appendix C of Ref. [18]. 

The zero neutrino (3(3 matrix element has, as mentioned above, contributions from the two components of the 
160 Dy wave function: 

M « = a M&(N„N V ) + b M^(A„N„) (35) 

As said before, these expressions are to be supplemented by other contributions, M^ l ,(N 7T A l ,) and M a !y (A 7r A„), when 
the large model space (of subsection III.B) is used to construct the initial and final wave functions. The expression of 
the additional terms is similar to the ones of the above equations and they are omitted for the sake of brevity. The 
matrix element of the /3/?oi/ mode is thus given by 

^E^^ M o,(k,l), (36) 

k,l 

where the indexes k and / denote the set of quantum numbers needed to specify the occupations and irreps included 
in the wave functions. 



VI. THE (3(3 OF GD 

In this section we study the two neutrino and zero neutrino modes of the double beta decay of 160 Gd to the ground 
state of 160 Dy. 

In the restricted configuration space of subsection III. A, the ground state of 160 Dy, defined in Eq.(8), is a linear 
combination of two states having different occupation numbers. The energy difference between these states, estimated 
from the difference in their deformed Nilsson single particle energies, is AE = 1.71 MeV. The pairing mixing between 
them, using the interaction strength = 21 /A MeV, is h pair — 0.865 MeV. With these matrix elements, the 
diagonalization of Eq. (9) yields the amplitudes a = 0.923, b = 0.385 of the wave function of the ground state of 
160 Dy. 

The (3(3 2v matrix element is suppressed by a factor 6, compared with the allowed decays. It implies that the 
(3f32 U half-life is an order of magnitude (1/& 2 ) larger than in other nuclei with similar Qpp values. The energy 
denominator takes the value £ = 12.19 MeV. The two neutrino (3(3 matrix element is M§J '( 160 Gd^ 160 Dy) = 0.0455 
MeV -1 . The phase space integral is G GT = 8.028 x 10~ 20 MeV 2 yr" 1 , using gA/gv = 1-0. The estimated (3(3 2v half-life 
is 

tH 2 ( 160 Gd -» 160 Dy) = 6.02 x 10 21 yr- (37) 

The contribution of the different angular momentum J to the matrix elements M^ V (N^ N„, J), in Eq. ( 35), are 
shown in Table III. 

It is remarkable that the J=0 channel largely dominates both the Fermi (F) and Gamow- Teller (GT) transitions. 
The 3 = 2 channel tends to reduce the transition matrix elements by about 20 to 30 %, For the Fermi matrix elements, 
the J = 4 channel also contributes noticeably. Fermi and Gamow- Teller matrix elements add coherently due to the 
sign inversion in Eq. (26). 

The transition matrix elements from the different neutron single particle angular momentum j v in the matrix 
elements MQ U (A 7T N u ,j u ), are presented in Table IV. In this case two neutrons are annihilated in the normal parity 
orbitals Pi/2,P3/2, /5/2, /V/2> ^9/2 and two protons are created in the intruder orbit hn/ 2 - 

As seen in Tabic IV, both for the Fermi and Gamow- Teller matrix elements the terms add coherently. While in 
the Gamow- Teller case the /ig/ 2 — » ^11/2 transition clearly dominates, in the Fermi case all transitions amplitudes are 
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2 
3 
4 
6 

Sum 



-0.11746 
0.03792 

-0.00131 
0.02090 
0.00312 

-0.05683 



GT 
0.25588 
-0.05016 
-0.00011 
0.00454 
0.00501 
0.21516 



M 0v (N n N v ) = 0.27199 
TABLE III: MSi,(N„N v , J) for the (3(3 0v of : 

iu F GT 



'Gd. 



Pl/2 
P3/2 
fc/2 

fr/2 
hg/2 
Sum 



0.01552 
0.02105 
0.02168 
0.06341 
0.04310 
0.16477 



-0.06792 
-0.05359 
-0.16133 
-0.11847 
-1.16598 
-1.56729 
-1.73206 



M 0v (A„N v 

TABLE IV: M^{A„N v ,j v ) for the (3(5 0v of 16() Gd. 



comparable. As in Table III, Fermi and Gamow- Teller final matrix elements add coherently. In absolute value, the 
allowed (3(3o v transition matrix element Mo v (A v N l ,) is a factor 6 larger than the forbidden one M 0v {N^N lJ ). The final 
(3(3q v matrix element is 

M „( 160 Gd -» 160 Dy) = aM Qv {N n N v ) + bM Qv {A K N v ) = 0.251 + 0.668 = 0.919. (38) 

The zero neutrino phase space integral for the (3(3q v of 160 Gd is Gq v — 1.480 x 10~ 14 yr -1 . Expressing the Majorana 
mass of the neutrino (m„) in units of eV , the calculated zero neutrino (3(3 half- life is 

r Y 2 ( 160 Gd -» 160 Dy) • (m„) 2 = 2.09 x 10 25 yr- (39) 

This half- life is a factor 20 larger than the one reported in [40]. The difference is easily understood, given the fact 
that the present model takes explicitly into account the nuclear deformation, while in [40] the spherical QRPA was 
used to obtain a crude estimations of the half-lives of all potential (3(3 emitters. 

As mentioned above, the parameter AE, which strongly influences the pairing mixing, is taken from the deformed 
Nilsson single particle energies. To estimate in which extent changes in this parameter affect the predicted doublc- 
beta-decay half-lives, we have allowed it to vary from 1.1 MeV to 2.5 MeV, covering most of the physically reasonable 
range. The results are listed in Table V. 

The mixing parameter b and the double-beta-decay matrix elements and half-lives are very smooth functions of 
AE. The two neutrino double-beta-decay half-life varies between 4 xlO 21 yr and 10 xlO 21 yr, around the predicted 
value of 6 xlO 21 yr. The calculated half-life of the neutrinolcss double-beta-decay is less dependent upon the mixing 
induced by pairing and it varies in the range 1.5 - 2.9 xlO 25 yr. 

Concerning the results obtained in the enlarged space of subsection III.B, we are presenting them in Table VI. The 
partial contributions to the (3(3i v and (3(3q v decay processes, for transitions between the components of the initial and 
final wave functions, are listed in this table. We are indicating, also, the character of each of the transitions, with 
regard to the type of orbital, abnormal (A) or normal (N) , of the neutrons and protons involved in the decay. 

For the case of the (3(3i v decay mode, there are four active configurations and each of them add coherently to the 
final matrix clement. They are partially suppressed by their amplitudes and the final matrix element is of the order of 
0.086 MeV _1 . This value is about twice the one obtained by using the small configuration space of subsection III. A. 
Consequently, the calculated half- life, which is of the order of 1.68 xlO 21 yrs, is shorter than the one obtained in the 
small space. Concerning the (3(3q u decay mode, the results shown in Table VI indicate that there is an interference 
between configurations where both nucleons are in intruder orbits and those where the proton occupies an intruder 
orbit. For this channel the resulting matrix element is of the order of 0.293, a value which is about a factor three 
smaller than the one obtained in the small configuration space. The predicted half-life 

To/ 2 ( 160 Gd -» 160 Dy) • (m„> 2 = 2.05 x 10 26 yr, (40) 
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AS[MeV] 


b 


M£ T [MeV _1 l 


r, 1/2 fl0 21 yrl 


M 0v 


Tn 10 vr 


1.10 


0.481 


0.0568 


3.86 


1.072 


1.53 


1.20 


0.464 


0.0547 


4.16 


1.044 


1.62 


1.30 


0.447 


0.0527 


4.48 


1.017 


1.70 


1.40 


0.431 


0.0508 


4.82 


0.992 


1.79 


1.50 


0.415 


0.0490 


5.18 


0.967 


1.89 


1.60 


0.401 


0.0473 


5.57 


0.943 


1.98 


1.70 


0.387 


0.0456 


5.98 


0.921 


2.08 


1.71 


0.385 


0.0455 


6.02 


0.919 


2.09 


1.80 


0.374 


0.0441 


6.41 


0.900 


2.18 


1.90 


0.361 


0.0426 


6.86 


0.879 


2.28 


2.00 


0.349 


0.0412 


7.34 


0.860 


2.38 


2.10 


0.338 


0.0399 


7.83 


0.841 


2.49 


2.20 


0.327 


0.0386 


8.36 


0.824 


2.60 


2.30 


0.317 


0.0374 


8.90 


0.807 


2.71 


2.40 


0.307 


0.0363 


9.47 


0.791 


2.82 


2.50 


0.298 


0.0352 


10.06 


0.776 


2.93 



TABLE V: The mixing parameter b, the double-beta-decay matrix elements and half-lives are listed as functions of the parameter 
AE. 



Transition 


channel 


2v>PI3 mode vfiP 


ii 




fi 


N n A v 


not allowed 


-0.231 


k 




h 


N^N V 


not allowed 


0.272 






h 


A^A^ 


not allowed 


-2.122 


ii 




h 


AkN„ 


0.118 


1.732 


ii 




h 


N^N V 


not allowed 


0.315 


ii 




h 




0.122 


1.787 


iz 




h 


At; Au 


not allowed 


-2.122 


h 




h 




0.118 


1.732 


i 4 




h 




0.122 


1.787 



TABLE VI: Matrix elements for the (3f32v and /3/3o v decay modes. Each arrow shows the participant configurations, of the initial 
and final wave functions, which are listed in Tables I and II. The transitions not included in the table are trivially forbidden 
because they imply the change of the states of more than two nucleons. The /3/32 V transitions listed as "not allowed ", are the 
ones which are forbidden by the selection rules discussed in section IV. The matrix elements of the (3/32v mode are given in 
units of MeV -1 . 

is an order of magnitude larger than the one obtained in the small configuration space. 

The above presented results can be summarized by noticing that the effect of including occupations others than the 
most probable one is less crucial for the two-neutrino mode than for the case of the neutrinoless double beta decay. 
Nevertheless, the predicted two-neutrino double beta decay of 160 Gd is still suppressed, as compared to other double 
beta decay emitters. In this respect, the results of the present calculations are an improvement of earlier ones [17], 
where claims about a suppression of the two-neutrino mode have been made. Here, we have used a larger configuration 
space, as explained before, instead of a single configuration. The two neutrino double beta decay in 160 Gd is hindered 
by nuclear structure effects, and the predicted half-life is of the order of 10 21 ( 22 ) yr, depending upon the model space. 
The zero neutrino double-beta-decay half-life is at least three to four orders of magnitude larger. In view of these 
predicted values, we are confident that the planned experiments using GSO crystals [24] would definitely be able to 
detect the flfcv decay of 160 Gd, and could establish competitive limits to the (3/3o„ decay. The background suppression 
due to a large flfcv half- life would be effective, although not as noticeably as was optimistically envisioned in [24]. 

Results about selection rules in other deformed double beta decay emitters are reported in [45]. 
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VII. CONCLUSIONS 

In the present paper we have studied the fifov and (3(3q v decay modes of 160 Gd to the ground states of 160 Dy. The 
transitions have been analyzed in the context of the pseudo SU(3) model. 

The energy spectrum and electromagnetic transitions in 160 Gd and 160 Dy have been studied in detail, in previous 
works, using the pseudo SU(3) model and a realistic Hamiltonian. Ground state wave functions were built as linear 
combinations of the pseudo SU(3) irreps associated with the larger quadrupole deformations, in a model space with 
fixed occupation numbers in normal and unique parity orbitals. Nucleons occupying intruder orbits were frozen. The 
pseudo SU(3) leading irrep typically carries 60% of the total wave function. 

In the present contribution the mixing of different occupation numbers in the 160 Dy ground state wave function 
was studied. Only leading irreps, for each occupation, were considered in the calculations. The mixing induced by 
the pairing interaction makes possible the two neutrino double beta decay of 160 Gd, which is forbidden when only the 
most probable occupation numbers are used. 

Explicit expressions are presented for the pairing mixing, and for the flfiiv and (3(3q v nuclear matrix elements in the 
present pseudo SU(3) approach. The estimated (3/3 half-lives are larger than those obtained using a spherical QRPA 
model, and the results suggest that the planned experiments would succeed in detecting the Pfov decay in 160 Gd, and 
in setting competitive limits for the zero neutrino mode. 
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